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Abstract 

Using Grad's method, we calculate the entropy production and derive a formula for the second- 
order shear viscosity coefficient in a one-dimensionally expanding particle system, which can also 
be considered out of chemical equilibrium. For a one-dimensional expansion of gluon matter with 
Bjorken boost invariance, the shear tensor and the shear viscosity to entropy density ratio rj/s 
are numerically calculated by an iterative and self-consistent prescription within the second-order 
Israel-Stewart hydrodynamics and by a microscopic parton cascade transport theory. Compared 
with r]/s obtained using the Navier-Stokes approximation, the present result is about 20% larger 
at a QCD coupling as ~ 0.3(with rj/s ^ 0.18) and is a factor of 2 — 3 larger at a small coupling 
Us ~ 0.01. We demonstrate an agreement between the viscous hydrodynamic calculations and the 
microscopic transport results on r//s, except when employing a small Us- On the other hand, we 
demonstrate that for such small a^, the gluon system is far from kinetic and chemical equilibrium, 
which indicates the break down of second-order hydrodynamics because of the strong noneqilibrium 
evolution. In addition, for large (0.3 — 0.6), the Israel-Stewart hydrodynamics formally breaks 
down at large momentum pt ~ 3 GeV but is still a reasonably good approximation. 

PACS numbers: 24.10.Lx, 24.10.Nz,12.38.Mh, 25.75.-q, 05.60.-k 



* el@th.physik.uni-frankfurt.de 
''' Azwinndini. Muronga@uct.ac.za 

xu@th.physik.uni-frankfurt.de 
^ Carsten.Greiner@th.physik.uni-frankfurt.de 



1 



I. INTRODUCTION 



Recent experimental measurements on the elliptic flow parameter V2 at the BNL Rel- 
ativistic Heavy Ion Collider (RHIC) show a strong collectivity of the deconfined 

quark-gluon matter. The matter produced was thus specified as a strongly coupled quark- 
gluon plasma (sQGP) ^, 0, or as a perfect fluid fzj]. Further attempts to determine how 



imperfect t 
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le sQGP re ally is have drawn attention to transport coefficients like viscosity 
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id] and to the derivation and solution of viscous hydrodynamics 
24| . which is still a mathematical challenge. 



Most current viscous hydrodynamic equations are based on second-order Israel-Stewart 



kinetic theory 



25j. They are solved numerically using the given viscosity coefficients and 



initial conditions as well as parton and hadron equation of state. In particular, the shear 
viscosity to entropy density ratio t]/s is determined by comparing the elliptic flow from the 
viscous hydrodynamical calculations with the data at RHIC, as has been done recently in 
Refs. 26|, |27l], where the value rj/s ~ 0.1 was obtained. On the other hand, even though the 
early partonic phase may be well described by ideal hydrodynamics [r] = 0), the hadronic 
afterburning 28|| has a larger dissipative effect, which may be enough to slow down the 
generation of the elliptic flow and bring its flnal value into agreement with the data. 

Dissipative phenomena can be alternatively described in transport calculations solving 



Boltzmann equations of matter constituents 



29 



30 



31 



32 
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34 



35 



36|. This approach is 



applicable for investigations of such phenomena as thermalization, kinetic decoupling, and 
dynamics of high-energy particles in systems far from equilibrium, i.e., in a regime where 
the second-order viscous hydrodjTiamics breaks down 37|. 

Recently, an on-shell parton cascade Boltzmann Approach of MultiParton Scatterings 



34, 



38|,l39|], elliptic flow t>2 



40 



41 



42|, 



(BAMPS) has been developed to study thermalization 

and the energy loss {43=] of gluons produced in Au+Au collisions at RHIC energy. Also the 
generation and evolution of viscous shock waves are surprisingly well realized in BAMPS 
calculations [441]. The shear viscosity of the gluon matter at RHIC has been estimated from 



BAMPS calculations 



41 



42| within the Navier-Stokes approximation 12]. The authors 



found that to produce large V2 comparable with the experimental data, the gluon matter 
should have a.n rj/s between 0.08 and 0.2 constrained by details of the hadronization and 
the kinetic freeze out. This is in line with the dissipative hydrodynamic approach 26). 
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Perturbative QCD (pQCD) gluon bremsstrahlung gg ^ ggg is responsible for the low rj/s 
ratio and for the generation of large elliptic flow. 



Beyond the Navier-Stokes approximation, which has been used in Refs.|lll. Il2[|. we derive 
a new microscopic formula for the shear viscosity coefficient from the kinetic theory using the 
second-order Grad's method. This is one of the goals in the present article. The derivation 
follows Ref. {451 and is generalized for a particle system out of chemical equilibrium. 

Another goal is to elaborate on the breakdown region of the second-order viscous hydrody- 
namics. To do this we investigate the time evolution of a gluon matter in a one-dimensional 



expansion with Bjorken boost invariance 46| by solving the Israel-Stewart hydrodynamic 



equations [37j as well as by performing similar BAMPS transport calculations for compari- 
son. We quantify the deviation of the gluon distribution function from kinetic equilibrium 
and show the region with large deviation, where the applicability of the Israel-Stewart hy- 
drodynamics is questionable. 

The article is organized as follows. In Sec. [TTl we introduce theoretical framework for 
deriving viscosity from the kinetic theory using second-order Grad's method. We consider a 
massless particle system, which undergoes a one-dimensional expansion with Bjorken boost 



12 



invariance. A comparison with the results obtained by the Navier-Stokes approximation 
is given in Sec. IIIII Using the formula derived in Sec. [TTl we calculate the shear viscosity 
to entropy density ratio rj/s oi gluon matter: in Sec. IIVI an iterative and self-consistent 
approach is introduced to calculate rj/s from the Israel- Stewart hydrodynamics, whereas the 
results from BAMPS calculations are presented in Sec. |Vl For both hydrodynamic and 
transport calculations, deviations from kinetic as well as chemical equilibrium are shown 
and analyzed. Conclusions are given in Sec. IVI[ 



II. SHEAR VISCOSITY COEFFICIENT FROM SECOND-ORDER KINETIC 
THEORY 

Relativistic causal dissipative hydrodynamic equations can be derived from the kinetic 
theory by ^applying Grad's method of moments 47|. A detailed derivation is reported in 



Refs. 



45 



48l ] and a prescription for calculating transport coefficients is also given there. In 



this section we will follow Ref. 45| to derive an expression for the shear viscosity coefficient 



7] when the considered system is out of chemical equilibrium. 
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The basic equation of relativistic kinetic theory is the Boltzmann equation 



p^d,fix,p) = C[f{x,p)] (1) 



for a one-particle phase-space distribution function f{x,p)) = i ' ^j^ ^3 -- C[f{x,p)] denotes 



the coUision term, which accounts for sl. 



The entropy four-current is defined by 45|, |49 1 



1 microscopic interaction processes among particles. 



s'' = -J J^^oP'fi^^P) Mf{x,p)) - 1] . (2) 
The entropy production is then given by 

d^s^ = - j dwp^'d^f{x,p)\nf{x,p) = -j dw C[f{x,p)]\n f{x,p) (3) 



with the short notation dw — Trr^fTTT- 

We now assume that the deviation of f{x,p) from the equilibrium distribution feq{x,p) 
is small: 

f{x, P) = feq{x, p) (1 + p)) (4) 

where (t){x,p) <^ 1 and 

feq{,X,p)=\e —. (5) 

\{x) and T(x) denote the local fugacity and temperature, respectively. u^{x) is the hydro- 
dynamic four-velocity of the medium. Equation ([5]) is the standard form for Boltzmann 
particles. The derivation below can be easily extended for Bose and Fermi particles. In 
addition, we will restrict the following discussions to the case of massless particles (e.g., 
gluons) . 

We expand (j){x,p) up to second order in momentum, that is, 

= e(x) - e^(x)p^ + e^,y{x)p^p'' , (6) 
where the momentum-independent coefficients can be expressed in terms of the dissipative 



currents 11, and n^'^ denoting bulk pressure, heat flux and shear tensor [45|, |48|: 



e^i^ = ^2(3u^M,. - A^^)n - + CoTT^^ (7) 

= AiUfJl - BaQf, (8) 
e = AoU (9) 



with the projector A'^'^ = g'^'' — u^u^ and symmetrization operation Uf^^q^) = \{u^qy + Uyq^). 
The metric used in this work is g^'^ = diag{l, —1,-1,-1). In general, the dissipative fluxes 
are deflned as projections of deviations of the energy-momentum tensor T^'^ and particle 



four-current from their equilibrium form 45|, |48|] : 



n = -^A^.ST^^-' (10) 
qf" = A'^UpST"" - ^A'idN'' (11) 



5T<^^-> = Qa^^A^ + ^A^A^ - ^A^^A^^^^ dT'^^ (12) 



with the deflnitions A^^ = / dwp^'f , Tf"" = J dw^fp^f and 6T^"' = T'^^ - T^^", dNi" = 

We use the following local matching conditions on the energy and particle densities: 

e = eeg = — ^ (13) 
n = neg = (14) 

with the deflnitions for the densities e = u^^T^^Uy and n = u^N^. The local temperature 
simply follows as T = e/3n. The fugacity is then calculated via X = n/ {^T^). One obtains 
immediately u^ST^^Uu = and Ufj,6N^ = 0. The bulk pressure 11 from Eq. ( fTOl) then 
becomes 

n ~ {g^u - u^u,)6T^'' = 6T;^ = (15) 

for massless particles, since the energy momentum tensor is traceless in this case. Thus, 
e = according to Eq. ([9]). 



In the following, we will consider a one-dimensional Bjorken expansion [46|]. This implies 
that in the local rest frame, the distribution function f{x,p) is symmetric when transforming 
p to —p. Thus in the local rest frame, T"* = and A^* = 0, where i = 1, 2, 3. The heat flux 
q'^ ( [TT]) vanishes in the local rest frame because 

= gj^UpSTP" - u^'UydTP" - -gi^dN" + -u^m^^A^^ = Up5T«' - -^A^^ = . (16) 

3 3 3 

We obtain then e^p^ ~ q^p^ = [see Eq. iQ]. 

For a one-dimensionally expanding system, Eq. ([6]) thus reduces to 

0(a;,p) =e;..(x)p>^ (17) 



Putting / = /eq(l + 0) into Eq. ([3]) and using the linearization 

ln(l + 0) ^0 = e^,(xKp" (18) 

we rewrite Eq.(I3]) as 

d.s'' = - j dwC[f{x,p)]\nU,{x,p) - j dwC[f{x,p)]e,,p^^p^ (19) 
Using the formula ([5]) for /eg in the first term of Eq. f|T9l) one has 

- j dwC[f{x,p)] ln/,g(x,p) = - j dwC[f{x,v)] (InA-nX/^) 
= -InA I dwC[f{x,p)]+u^ ! dwp^C[f{x,p)]/T 



= -InA j dwC[f{x,p)] = - \n\df,Nr (20) 

For the second-last identity in Eq. (!20|) . we used the energy-momentum conservation: 
f dwp'^ C[f{x,p)] = df^ f dwp^p^f = dfjT^^ = 0. Equation ( l20l) describes entropy pro- 
duction due to particle production {d^N^ > for A < 1) and absorption {d^j^N^ < for 
A > 1). 

With the definitions 

P"'' = j dwp>'p''C[f{x,p)] (21) 

C = J dwC[fix,p)]=d^N'^, (22) 

which are the 2nd and the 0th moment of the collision term the entropy production in 
Eq. (fT9!) can be now written in a more compact form 

a^s^ = -ClnA-e^,P'^^ (23) 

In general, the entropy production in an imperfect fluid can be expressed by the positive 



definite form 



25 



50 



5l| 



d^s^ = -Jin A + {CT)-^U' - {^T)-\^q" + (2r/T)-V,^7r"^ , (24) 



where (, k, and rj are non-negative coefficients denoting the bulk viscosity, heat conductivity 



and shear viscosity, respectively. J = d^N'^ is the source of particle production [50|, l51| and 
is identical with C (1221) . For a chemically equilibrated system J vanishes. Comparing Eq. 
([23D to ([21D we find 

-e^,P'^^ = (277T)-V„;37^"^ (25) 



because in our case 11 = and qaq" = as discussed above. The expression (125!) is exactly 



the same as obtained in 



45j and describes entropy production due to shear viscous effects. 



We then obtain the final expression for the shear viscosity coefficient 



V 



a/3 



al3 



(26) 



The last identity is due to the fact that vanishes in the local rest frame and thus 
U{^qu)P^^ = 0. We note that the derived formula (12^ is an approximate expression of 
the true shear viscosity. We call the "second-order" shear viscosity, because we have used 
terms up to second order in momentum for (j){x,p) [see Eq. (Q]. 

To calculate Cq we go to the local rest frame, i.e., u^^ = (1, 0, 0, 0), where 



TT 



ST^" = T^" - T^^ = ea/3 j dwp^p''pyU,{x,p) 



(27) 



is valid according to Eqs. f|T2|) and (fT71) for a (0+1) dimensional expansion. In this frame 
e^p [see Eq. ([7])] reduces to 

ea/3 = CQ-KaP ■ (28) 

Calculating the integrals in Eq. fl27|) with f^q = Xe~^^^ gives 



;i-Co40ATV7r2)7rO^=0, j = 1,2,3 
;i - Co8ArV7r2)7r*-'' = , i,j = 1,2,3. 



(29) 
(30) 



We have used the fact that vr^'^ is traceless and 7r°° = due to the matching condition (fT3|) 
and = e in the local rest frame. For a system undergoing a one- dimensional Bjorken 
expansion, i.e., in a (0+1) dimensional case, all off-diagonal elements of T'^" - and thus tt^" 
as well - vanish in the local rest frame, particularly T^^ = vr"-' = 0,j = 1,2,3. Thus we 
obtain 



Co 



TT 



(31) 



8AT6 ■ 

If the third spatial coordinate is chosen as the expansion axis, we have T^^ = T^^, and in 
the local rest frame the shear tensor takes the form 



TT' 



/o o\ 



-f 



-- 

2 ^ 




\0 tt/ 



(32) 



which is also given in 51i- We thus obtain 



vr^.vr^'^ = (33) 



e^.P'" = Con.^P'^'' = ^ {3P'' - P°0) , (34) 



where we have used P^^ + P^^ = P^^ — P^^, because P^*^ is traceless following from the 
definition fl2Tl) . Putting Eqs. fl27j) and flMl) into fl26l) gives the shear viscosity coefficient for 
a (0+1) dimensionally expanding system of massless particles: 



37f . -T 



TT 



^ 2TCo(3P33 - poo) ^'^pss.ipoo- (^^) 

For the last identity, we have used the matching conditions (HM and Eq. (pTi) . 

The energy density e, the temperature T and the shear component vr in a (0+1) dimen- 
sional expansion can be calculated by solving viscous hydrodynamic equations with a given 
value of shear viscosity t]. If vf is known, the distribution function / is known too [see Eqs. 
(jl]), (|T7j) and fl28|) ]. One can thus evaluate P"" and P^^ according to their definitions fl2Tl) . 
Then r/ can be calculated employing Eq. ( !35l) . In sSec. IIVI we will introduce an iterative 
and self-consistent prescription to calculate the second-order shear viscosity. 

On the other hand, / can be obtained by solving the Boltzmann equation ([T]) directly 
employing transport simulations. Then rj can be easily extracted using Eq. (!35|) . Such 
calculations will be presented in section |Vl The results will be compared with those obtained 
in Sec. IIVI As it turns out, a ratio of mean transport free path to expansion time being 
larger than unity and the variance of (f){x,p) being larger than unity will possibly indicate 
the breakdown of the second-order viscous hydrodynamics. In this regime the validity of 
fl55]) is also questionable. 



III. COMPARISON TO SHEAR VISCOSITY FROM NAVIER-STOKES AP- 
PROXIMATION 



12|, the shear viscosity coefficient was derived assuming the Navier-Stokes ap- 



In Ref. 

proximation 

7c^''^ = 2r]V<^u''> . (36) 

It reads 

^1 {E/3-pI/E)) 1 
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where 

otr _ Idw^C[f]-{pl/E^)JdwC[f] 
^ n{\-{vVE^)) ^'^^ 

is the total transport colhsion rate, which was introduced in 38|. All integrals are expressed 
in the local rest frame. () denotes the average over particle momentum. 

Equation (126!) can be used to calculate the shear viscosity if the shear tensor t^^^ obeys 
the Israel-Stewart equation [51] 



(39) 



where a^^ = V^^u^^ and denotes the relaxation time [see also Eq. fjTr|) below]. Equation 
(!39|) is more general than (136|) in the first-order (Navier-Stokes) theory. 
If we define 

p33 _ IpOO 

then the shear viscosity from the Grad's method ( !35l) can be rewritten to 

t\e/?,-pIie) 1 

E^S™/ ^ ^ 

where we have used vf = T^^ - T^s ^ ^"33 _ 1^°° = n{pl/ E - E/3). Remember that Pf"" is 
the second moment of the collision term [see Eq. fl2T]) ]. The expression (j^Tl) is similar to Eq. 
( !37l) except for the term |(9t(lnA), which indicates that chemical equilibration contributes 
explicitly to the shear viscosity in the Navier-Stokes approximation rather than in the Israel- 
Stewart approach. 

In the next section, we calculate the shear viscosity in a gluon system within the Israel- 
Stewart approach and compare the result with that obtained using the Navier-Stokes ap- 
proximation [12]. 



IV. CALCULATION OF SHEAR VISCOSITY IN A GLUON SYSTEM: AN ITER- 
ATIVE AND SELF-CONSISTENT PRESCRIPTION 

In this section we want to calculate the shear viscosity to the entropy density ratio rj/s 
for a gluonic system, which undergoes a one-dimensional expansion with Bjorken boost 
invariance, i.e., a (0+1) dimensional expansion. 
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A. Prescription 



For a (0+1) dimensional case the shear tensor vr^j, in the local rest frame is given by Eq. 
Then the gluon distribution function in the local rest frame reads 

f{x,p) = Ae-f [1 - C.TTipl - p'j2)] (42) 

according to Eqs. (jl]), (fT7|) . ( l28l) and ( l32l) . If 7f, T and A are known, the shear viscosity 
T] can be calculated according to Eq. (^5^, where P^'^ are evaluated by Eq. (pTi) via Eq. 
f H2]) . Note that for the case of a gluonic system the value of rj has to be amplified by the 
degeneracy factor of gluons dc = 16. We thus define rjg = dcr]. In addition, the gluon 
entropy density is given by 



Sg U^S^ 



-de J dwpofix,p)i\nfix,p)-l)^iA-lnX)ng-^^, (43) 



where Ug = (ic-AT^/Tr^ and vf^ = dcif are the gluon number density and the gluon shear 
component, and we have used the approximation ln(l+0) ^ for small = —Cqtx{p1—p1/2). 
We note that can be larger than unity for large momenta. In these cases, the expansion 
[also for Eq. (fT9l) ] fails. On the other hand, the distribution function f{x,p) becomes very 
small at large momenta. The effect of the invalid expansion on the integrated quantity Sg is 
thus negligible at this point. 

In principle, = (/ — /eg) / /eg gives the relative deviation from kinetic equilibrium. How- 
ever, is also a function of momentum. The average (0(x,p))eg over momentum distributed 
in equilibrium, i.e, using f{x,p) in zeroth order of vf, is obviously zero. We introduce the 
variance cr^ = a/ (0^)eg as the quantity determining the deviation from kinetic equilibrium 
and we find 

(44) 



^5 1 



4 eg 



where Cg = SugT is the gluon energy density. 

If the deviation from the local kinetic equilibrium is sufficiently small, then the dynamical 
expansion in a (0+1) dimensional case can be well described by the Israel-Stewart (IS) 
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viscous hydrodynamic equations |25|, l37|, |45|, |48|, l5l|, l52| : 



dUg 


_ng 


dr 


T 


deg_ ^ 


4 eg 


dr 


3 r 


dTig 


_^ _ 


dr 





r 
1_ 

2 ^ 



p2 dr^T' 



2 1 



3/32r 



(45) 
(46) 
(47) 



where P2 = 9/(4ec,) and = 2P2Vg denotes the relaxation time. Equation ( H7I) is just Eq. 
(!39|) expressed in the local rest frame using the hydrodynamic velocity u'^ = -(t,0,0, z), 



where r = ^/t"^ — z^. In derivation of Eq. fl39|) . which is discussed in Ref. [51,], only terms 
of second order in gradients and dissipative flux vr^,^ have been included. If in Eq. (jH]) 
is larger than unity, further terms containing cr| ~ i^/^Y ~ (/527r)^ are no longer small 
enough anymore to be ommited in derivation of Eq. ( l39i) and thus in Eq. ( 1471) as well, i.e., 
a higher order hydrodynamic equation is needed. Thus the value of cr,^ is an indicator for a 
breakdown of second-order hydrodynamic theory. 

Equation (H5|) for the gluon density can be easily solved: 



rigir) = ng{To) — 



(48) 



which is identical with the result from ideal hydrodynamics. On the other hand, the energy 
density decreases slower than in ideal hydrodynamics due to the viscous effects: 

'ro\« . . 4 



(49) 



Thus we obtain the gluon fugacity 

_ ngir) 



A(. 



n 



eq. 



Ao(^)'"^'<Ao, 



(50) 



where Aq = A(ro). The system will be continuously out of chemical equilibrium during 
the expansion, even if it is initially at local thermal equilibrium (Aq = !)• The larger the 
viscosity, the smaller is the value of ^ and the faster is the decrease of the fugacity. Inclusion 
of production and annihilation processes such as the gluon bremsstrahlung and its back 
reaction {gg ^ ggg) makes chemical equilibration possible and thus, of course, Eq. fHS]) has 
to be modified! However, in this work we will use Eq. fl45p without any modifications. The 
derivation of new and altered equations and their solutions will be given in a forthcoming 
publication [53|. 
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One can solve Eqs. (H6i) and (H7I) . if the initial values of n^, eg, 7tg and also the value of 
the shear viscosity rjg are given. On the other hand, to calculate rjg using Eq. fl35|) via Eq. 

we need n^, e^, and vfg. It is obvious that an iterative algorithm has to be developed 
to calculate Ug, Cg, fCg and r]g self-consistently. This algorithm is as follows: 

1. We solve Eqs. fj45l) - P7j) with a guessed value of rjg. The guessed value can be chosen 
arbitrarily because the final result does not depend on it. rig/ug is assumed to be a 
constant of time. 

2. The obtained ng^r), Cgir) and vrg(r) at a time r are used to calculate ?7g(r) according 
to fl35|) . We calculate first the moments P^^ using f{x,p) in Eq. (H2|) with given 
ng(r),eg(r) and vfg(r). 

3. We turn back to step 1. The value of rjg^r) is used to solve Eqs. (H6|) and (H7|) again. 

Iterations will continue, until the relative deviation of rig from the previous one is sufficient 
small. The iterative procedure allows to calculate vr(r), e(r) and n(r) as well as ?7/s(r) in 
a consistent way for given interactions. We note that if Tjg/ng is strongly time dependent, 
further iterations will be required to account for this time dependence. A refined algorithm 
will be presented in 

To obtain rjg, P^^ has to be first evaluated by fl^Tl) via fH21) . P^^ is a second moment of 
the collision term and thus is determined by gluon interactions considered. The compact 



forms of the collision terms can be found in 



34| . In this article elastic {gg gg) as well 



as bremsstrahlung {gg ^ ggg) processes inspired within perturbative QCD are responsible 
for the gluon dynamics. The differential cross section and the effective matrix element are 



taken as in Refs. 



34 



39|: 



'-^—^1' = — (qi+m^)2 ki[(k^-q^)'2+m|,] ^^^^^^ ' ^'2) 
where g^ = ATrag- The Debye screening mass 

ml, = dcTras J dw Ncf{x,p) (53) 

with A^c = 3 is applied to regularize infrared divergences. Although gg ^ ggg processes are 
considered, they contribute only to the shear viscosity but not to chemical equilibration, 
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because as mentioned above, particle number conservation is assumed at present to derive 
Eq. fH5l) ). Improvements will be done in a forthcoming publication 53|. 



B. Results 

Figure [U^a) shows rjg/ Sg as a function of the expansion time for two values of the coupling 
constant = 0.05 and 0.3. The initial gluon system at tq = 0.4 fm/c is assumed to be in 
thermal equilibrium with a temperature of Tq = 500 MeV. Each of the results indicated by 
the symbols in Fig. [T]is obtained by about 40 iterations with a guessed value of ?7g(guessed) = 
0.5 s^''. From Fig. [I](a) we see that the ratio fjg/sg is almost constant in time for as = 0.3, 
whereas for as = 0.05, f]g/sg increases moderately. The assumption underlying the iterative 
algorithm that Tjg/ng ^ Arjg/sg does not depend on time is justified accordingly. One finds 
that rjg/sg ~ 0.18 for a coupling of = 0.3 and rjg/sg ~ 3 for = 0.05. 

The results for the gluon fugacity (obtained from the solution of Eqs.( l45i) -( l47l) ) depicted 
in Fig. [H^b) show a strong time dependence. The smaller the value of as, i.e., the larger 
the fjg/sg, the faster is the deviation from the chemical equilibrium. This quantitatively 
demonstrates the consideration from above [see Eqs. fl50l) ]. 

When putting Eq. (142|) into Eq. (pTj) one realizes that P^'^ ~ X^CqTT ~ Att in leading order 
of 71. Thus rjg does not depend on A. Secondly, from Eq. fH5]) we obtain Sg/T^ ~ A(l — In A). 
Thus, Tjg/sg ~ 1/A(1 — In A) and will increase slower than a logarithmical behavior when A 
decreases: a stronger decrease of A (comparing the result for = 0.05 with that for = 0.3 
in the lower panel of Fig. [1]) will lead to stronger increase of Tjg/sg, as seen in the numerical 
results shown in Fig. [T](a). 

Figure [2t^a) shows the deviation from kinetic equilibrium, from Eq. (jS]), as a function 
of time scaled with the initial time. For as = 0.3 the value of starts at zero (equilibrium), 
increases until Stq and then relaxes to zero. The system first evolves out of equilibrium and 
then relaxes back to equilibrium. On the contrary, increases continuously when employing 
a much weaker (and unphysically low) coupling as = 0.05. In this case the system is always 
out of equilibrium. To explain the different behaviors we define Rqe as the ratio of the 
mean transport free path, l/X]-^Grad defined by Eq. ( l40l) . to the Hubble-like expansion 
time scale r: 



Roe = — 



1 



(54) 



Grad 
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Our concept of Rqe is similar to that introduced in [37|, where the authors demonstrate that 
the ratio of expansion time to the mean free path controls the deviation from equilibrium. 
For a fixed Tjg/sg the mean transport path A*'' = l/Xl-^crad changes with time. At full 
equilibrium A**^' ~ 1/T ~ r^/^ and thus A*^'/r ~ r^^/^. If Rqe is larger than unity, the 
system starts to depart from equilibrium. If Rqe is smaller than unity, the system relaxes 
to equilibrium. Roe{t) is shown in Fig. [2]^b). With as = 0.05 the system evolves far away 
from equilibrium and the evolution is dominated by free streaming. The ratio Rqe is a 
measure of the ability of the system to relax to kinetic equilibrium. For = 0.05 kinetic 
equilibration is not possible for the timescales shown. The regime for which the system can 
not come close to kinetic equilibrium is for the coupling = 0.1 — 0.2 corresponding to a 
shear viscosity to entropy density ratio rj/s = 0.8 — 0.4. 

In addition, is larger than unity at r > 3ro for as = 0.05. The true entropy density 
Sg should be smaller than that estimated according to Eq. (H3|l . because the expansion 
ln(l + 0) ^ is not valid any more for large 0. The derivation of the shear viscosity in Eq. 
(!35l) becomes questionable as well, since the same expansion is used to obtain the entropy 
production ( |T9|) . 

Finally, in Figs. [3] and H] we compare the results on rjg/sg from the second-order (IS) 
kinetic theory with those presented in Ref. [12|] using the Navier-Stokes approximation. 
The sohd (dotted) curve in Fig. [3] depicts the contribution of gg gg {gg ^ ggg) to rjg/ Sg 



obtained in [l^. The solid (dotted) curve with symbols depicts the results from the present 
calculations at r = 2ro, at which the system is still near thermal equilibrium. We see that 
the results following from the second-order expansion are mostly larger than those based on 
the Navier-Stokes scheme, both for gg gg and for gg ^ ggg processes. At (unphysical) 
small as the difference between the results is given by a factor of 2 — 3. In particular, 
the difference between the second-order and the Navier-Stokes results for bremsstrahlung 
99 ^ 999 is bigger than that for elastic gg gg process. At large ag the gg ^ ggg 
processes play a dominant role (compared with gg gg) in lowering rjg/sg, whereas at small 



as this dominance becomes weaker [ll|. In Fig. H] the results onrf/s implementing both 



elastic and inelastic processes are shown for the physical region of as- Here the difference 
between second-order and Navier-Stokes based calculations is approximately 50%(q;s = 0.2)- 





20%(a, = 0.3)-0%(a, = 0.6). 
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V. CALCULATION OF SHEAR VISCOSITY IN A GLUON SYSTEM: TRANS- 
PORT SIMULATIONS EMPLOYING BAMPS 



In this section, we solve the Boltzmann equation for gluons using the parton cascade 
Boltzmann Approach of MultiParton Scatterings and repeat the task in the previous sec- 
tion to calculate the shear viscosity to entropy density ratio Tjg/sg in a Bjorken-type one- 
dimensional (0+1) expansion. We calculate rig and Sg according to Eqs. (!35|) and (143|) by 
extracting P^^ , fCg, Ug and Cg from the transport simulations. 

The partonic cascade BAMPS which was introduced in 3^, ISSj has been applied for a 
(0+1) dimensional expansion to study thermalization of a color glass condensate potentially 

I I 

produced in ultrarelativistic heavy ion collisions [39]. We take the same numerical setup 



for BAMPS as considered in 



39| . The initial condition and interactions of gluons are the 



same as given in the previous section. In the parton cascade calculations, different from 
calculations using the viscous hydrodynamic equations (H5|) -f l47|) . the inelastic gg ^ ggg 
processes lead to a net particle production or absorption, i.e., d^N^^ ^ 0, which drives the 
chemical equilibration. 

We note that particle number changing processes are implemented in BAMPS, whereas 
the particle number was considered to be constant in previous section. Therefore we are not 
able to make a direct comparison between BAPMS results and those calculated by solving 
Israel-Stewart equations. 

Figure [5] shows rigjsg extracted within the space time rapidity interval 77^ G [—0.1 : 0.1], 
where ?7s = | ln[(t + 2;)/(t — 2;)]. When comparing these results with those shown in the upper 
panel of Fig. [T]we find that they are almost the same for = 0.3, whereas for = 0.05 the 
increase of Tjg/sg is slightly weaker in BAMPS calculations than in viscous hydrodynamic 
ones. The reason for this difference is the different behavior of the gluon fugacity (remember 
that f]g/sg ~ 1/A(1 — In A)). The gluon fugacity extracted from BAMPS is shown in Fig. [61 
Its value is larger than that shown in the lower panel of Fig. [H because ongoing chemical 
equilibration is realized in the BAMPS calculations. 

The kinetic equilibration is demonstrated in Fig. El^a) via the variance o"^ and in Fig. 

2 

[8] via the momentum isotropy Q{t) =< ^ >. The results on are similar to those in 
Fig. [2] and can be well understood by the out-of-equilibrium ratio Rqe shown in Fig. [7](b). 
For as = 0.3 the transport mean free path is shorter than the expansion rate whereas for 
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as = 0.05 the evolution is dominated by expansion. Momentum isotropization, shown in Fig. 
m is practically restored for as = 0.3 at later times, whereas for as = 0.05 this restoration 
is not possible. Here again the differences between transport and viscous hydrodynamic 
calculations stem from the different time evolution of the gluon fugacity. To make fair 
comparisons, modifications in the hydrodynamic equations will be done in the near future 



531] to take into account the chemical equilibration. 

As pointed out already in the previous section, the parameters for which the system 
cannot come close to kinetic equilibrium are the couplings of ~ 0.1 — 0.2 corresponding 
to a ratio 77/5 ~ 0.8 — 0.4. For such parameters, the ratio X^^/t is of the order of 1 at late 
times and the system becomes highly diffusive and viscous. 

Finally, in Fig. we investigate deviations from equilibrium of the gluon distribution 
in BAMPS calculations at large momentum. Figure [9] shows the non-equilibrium part of 
the transverse spectrum (normalized to the equilibrium spectrum) dN^J'(j)t%>T<ir]) ~ ^ from 
BAMPS calculations and the quantity < >p^ (7r,T, A) = / feq(l){7c,T, \)dp^/ J feqdpz [with 
(f){x,p) = oj^T^ipl — |Pt)], which is the analytically calculated second-order contribution 
to the transverse spectrum, as a function of the transverse momentum pt aX t = Atq. 
The average < (f) >p^ is calculated using vr, T, A extracted from the particular BAMPS 
calculations. The comparison of ^'^/(pj'^pj^^)^ _ i gj^^ < from Fig. [9] shows that for 
as = 0.05 the distribution function in BAMPS contains contributions higher order in p^ 
and TT and thus the second-order ansatz (fTTl) is not sufficient to describe the evolution in 
BAMPS. In contrast, for as = 0.3 the distribution function is reasonably good approximated 
by second-order kinetic theory over the shown momentum range. Thus we argue that the 
dependence of on 7f/(AT^) is stronger than given by ansatz (fT7|) . since 7f/(AT'^) ~ 
quantifies the strength of dissipative effects, which are stronger at as = 0.05. Inclusion of 
additional terms in Eq. (fT7j) would lead to a modification of the evolution equation for vf, 
which follows from the conservation law for the energy momentum tensor: = d^T^'^ ~ 
(9^/ p^p^ feqil + (p). If employing as = 0.05{ri/s ~ 3) for large pt > 2.3 GeV the variance 
< >p^ becomes larger than 1. For as = 0.3 this happens at pt > 2.75 GeV. For transverse 
momenta larger than these critical values the expansion /n(l + 0) ^ done to obtain Eq. ffTO]) 
is invalidated. Thus in the calculation of the entropy density (and entropy production) the 
ln(l -|- 0) term should be approximated by ln(l -(-0) ^ — ^&{pt — PTcrit), i-e. higher order 
terms should be taken into account in the integration over the momentum for p^ > prcru- 
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However, with rj/s = 0.18 this correction is less than 0.5%, which is due to the smallness of 
a^f). With rj/s = 3 the correction is ~ 6%. Thus for physical values of ry/s ~ 0.2 second-order 
hydrodynamics is valid, even though formally breaking down at large pt- In the unphysical 
regime ~ 3 higher order corrections are not negligible. This deserves future investigation 



53|. 



VI. CONCLUSIONS 

We have derived the shear viscosity coefficient rj from kinetic theory for massless par- 
ticle system undergoing a one-dimensional expansion with Bjorken boost-invariance. The 
derivation makes use of Grad's moment method 



45 



47( 1 and is based on an expansion of 



the distribution function up to second order in momentum. The final expression obtained 



in the present work is similar to the one based on the Navier-Stokes theory 12(], but the 



transport rate has to be calculated in a different way. How close the result obtained u sing 



Grad's method approximates the true value determined using the Kubo-Green formula 54 1 
will be studied and reported in a forthcoming publication. The values needed to calculate 
the shear viscosity [Eq. fl35l) ] are shear tensor vr^j^, the particle and energy densities e and n, 
the fugacity A and finally the second moments P^,y of the collision term from the underlying 
kinetic process. They can be calculated either using transport setup solving the kinetic 
theory or from dissipative hydrodynamic (Israel-Stewart) equations (H5l) -( H7I) . However, the 
IS equations themselves need the value of shear viscosity as a parameter. Thus we introduce 
a new iterative method that allows us to solve Israel-Stewart equations and calculate 7]/s as 
a function of time and coupling constant a.^- The results on rj/s calculated in the partonic 
cascade BAMPS and from IS theory are in a good agreement for physical coupling = 0.3. 
In this regime we obtain t]/s = 0.18. As a further demonstration even for unphysical small 
coupling as = 0.05 the difference between BAMPS and second-order hydrodynamic calcula- 
tions of rj/s is small. We obtain r^/s ~ 3 in this regime. At such small coupling rj/s increases 
slightly in BAMPS and somewhat stronger in hydrodynamic calculations. This increase is 
due to the intrinsic fugacity, which evolves differently in both calculations. 

Using IS theory, we calculate r]/s ratio for a system close to equilibrium as a function 
of as- For physical coupling ~ 0.3 the second-order result is approximately 20% higher 
than in calculations based on first order Navier-Stokes theory [l^. For as = 0.6 Tj/s = 0.08 
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within the Israel-Stewart and Navier-Stokes prescription. 

Deviations of hydrodynamic evolution from equilibrium are quantified in the present 
work introducing the variance of the nonequilibrium part of the distribution function. 
We demonstrate that its value is smaller than unity and later decreases with time at a 
physical coupling = 0.3 and thus our expression of rj is valid in this case. Here again 
hydrodynamic and BAMPS results are in good agreement. For small coupling as — 0.05 
hydrodynamics does not relax back to equilibrium and Grad's method becomes invalid. In 
BAMPS in this regime the deviation of from equilibrium is smaller, which is an effect of 
the ongoing chemical equilibration. The ability of the system to restore kinetic equlibrium is 
quantified by the ratio of the mean transport free path to the expansion time. We conclude 
that the second-order dissipative hydrodynamics is applicable in the regime ri/s ^ 0.2 which 
corresponds to values of ^ 0.3. At high momenta pt > 3 GeV it fromally breaks down, 
however for 77/s ~ 0.2 — 0.4 it is apphcable even for differential observables. For really high 
77/s ~ 3 the applicability of hydrodynamics certainly breaks down. For the intermediate 
regime 0.3 < rj/s < 0.8 a more detailed analysis in the comparison of microscopic transport 
description to dissipative second- (or even higher) order hydrodynamics is required. 

To make consistent comparisons between the kinetic transport model BAMPS and IS 
solutions we have to modify the hydrodynamic equation to take into account particle pro- 
duction and absorption. These calculations will be reported in a forthcoming publication. 
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FIG. 1: (Color online) (a) ry/s ratio and (b) fugacity A calculated by the iterative procedure 
described in the text for Og = 0.05 and Og = 0.3 at ten different time points, with initial time 
To = 0.4 fm/c, T(ro) = 500 MeV. The initial input value of r]/s is 0.5. 
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FIG. 2: (Color online) (a)Variance a^j, and (b) ratio Rqe calculated by the iterative procedure. 
To = 0.4 fm/c. 
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FIG. 3: (Color online) Ratio rj/s (contributions due to elastic and inelastic processes) ^as function 
of the coupling constant as- The result (solid line) is compared with results of Ref. (dotted 
line) 
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FIG. 4: (Color online) Ratio r]/ s 
is compared with result of Ref. 



(all processes) as function of the coupling constant . The result 
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FIG. 5: (Color online) Ratio r]/ s from the microscopic BAMPS simulation. Results are calculated 
using Eq. (f35]l for simulations with two different values for a^. tq = 0.4 fm/c 




FIG. 6: (Color online) Fugacity A = n/rieq from BAMPS calculation. Results are shown for 
simulations with different (constant) values of as- tq = 0.4 fm/c. 
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FIG. 7: (Color online) (a) Variance a^j, and (b) ratio Rqe calculated by BAMPS. tq = 0.4 fm/c. 
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FIG. 8: (Color online) Momentum isotropy < pi/E^ > calculated by BAMPS. tq = 0.4 fm/c 
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FIG. 9: (Color online) Nonequilibrium part of the transverse spectrum (normalized to the equilib- 
rium spectrum) ^'^/(Pr^Pjrf^)^ _ ^ from BAMPS calculations (lines with points) and the second- 
order contribution to the transverse spectrum / J'-^)'^?^ (lines) as function of at r = Atq 

J Jeqapz 

with 7f,T, A extracted from BAMPS. 
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